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APPENDIX A 

#include "main.h" 

#define XDIM 5121. 
#def ine XDIMR 5X2 
#define YDIM 4 80L 
#define BITS 6 
#define RMS VAL. 5.0 
#define NUM^NOISY 16 
#define NUM_DEMOS 3 
#define GRAD_THRE SHOLD 10 

struct char__buf { 

char filename [80] ; 

FILE *fp; 

f pOS_t £pOS ; 

char buf [XDIMR] ; 

struct uchar_buf { 

char filename [80] ; 

FILE *fp; 

f pOs_t fpos ; 

unsigned char buf [XDIMR] ; 

struct int_buf { 

char filename [8 0] ; 

FILE *fp; 

f pos_t fpos ; 

int buf [XDIMR] ; 

struct cortex_a { 

char filename [80] ; 
FILE *fp; 

f pOS_t fpOS ,- 

unsigned char buf [XDIMR] ; 



struct uchar_huf test_image; 
struct char_ buf snov^composite; 
struct uchar^buf distributed__image; 
struct uchar*buf terap_ image ; 
struct int_bu£ terap^wordbuf fer; 
struct int__buf tetnp_wordbuf f er2 ; 
struct uchar_buf snow^ images ; 
struct cortex_s cortex; 

int demo=0,- /* which demo is being performed, see notes */ 

int our_code; /* id value embedded onto image */ 
int found_code»0; /* holder for found code*/ 



int waitvbb(void) { 

while C (_inp (PORT_BASE) &8 ) ) ; 
while C \ (_inp<PORT_BASE) &8) ) ; 
return (1) ; 

} 

int grabb(vo±d){ 
waitvbb { ) ; 

_outp <PORT_BASE+l , 0) ; 
_outp (PORT^BASE, 8 ) ; 
waitvbb () ; 
waitvbb ( > ; 

_outp (PORT_BASE, 0x10) ; 
return (1) ; 
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} 

int livee (void) { 

_outp (PORT__BASE, 0x00) ; 
return(l) ; 



int live_video (void) { 

livee {) ; 
^ return (l) ; 

int f reeze_f rame (void) { 
grabb() ; 
return ( 1) ; 

} 

int grab_f rame (struct uchar_buf * image ) { 
long i; 

grabb() ; 

f setpos (image ->fp, &image->£pos ); 
f 3 etpo3 (cortex. fp, fccortex.fpos ); 
for (i«0;i<YDIM;i++) { 

f read (cortex . buf , ai2eof (unsigned char) , XDIMR, cortex, fp) ; 
^ f write (cortex. buf , sizeof (unsigned char) , XDXMR, image- >fp) ; 

livee () ; 
return ( 1) ; 

} 

int wait_vertical_blan3cs (int number) { 
long i; 

for { isO ; i <number ; i++ ) waitvbb ( ) ; 
re turn ( 1 ) ; 

} 



int clear_char_image (struct char_buf *charbuffer) { 
long i , j ; 
char *pchar; 
f pos_t trap_f pos ; 

f setpos (charbuf fer- >fp, acharbuf f er->f pos ); 
for (i-0;i<YDIM;i++) { 

f getpos (charbuf f er->fp, &trap_f pos ) ; 

pchar «=» charbuf fer- >buf ? 

f read (charbuf fer->buf, sizeof (char) , XDIMR, charbuf fer->fp ) ; 
for (j = 0/ j<XDIM; j++) *(pchar++) - 0 ,- 
f setpos (charbuf fer->fp, &tmp_fpos ) ; 

f write (charbuf fer- >buf f sizeof (char) , XDIMR, charbuf fer- >fp) ; 
return (1) ; 



int display_uchar (struct uchar_buf *image,int stretch) { 
unsigned char *pimage; 
unsigned char highest *= 0 ; 
unsigned char lowest = 255; 
long i, j ; 

double dtemp, scale, dlowest ; 
f pos_t tmp_f pos i 

if (stretch) { 

f setpos (image- >fp, &image~>fpos ); 

f read (image- >buf, sizeof (unsigned char) , XDIMR, image ->fp) ; 

f read (image->buf, sizeof (unsigned char) , XDIMR , image- >fp) ; 
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for (i~2;i< (YDIM-2) ;i+ + ) { 

f read < image ->buf , sizeof (unsigned char) ,XDIMR, image- >fp) ; 
p image - tiroage ->buf [3] ; 
for (j=3; j< (XDrM-3) ; j +-0 { 

if ( *pimage > highest ) highest « *piraage; 

if ( *pimage < lowest ) lowest = *p image; 

pimage++ ; 

} } 

if (highest « lowest ) { 

printf C something wrong in contrast stretch, zero 

contrast " ) ; 

exit(l); 

scale « 255.0 / < (double) highest - (double) lowest >; 

dloweat = (double) lowest ; 

f setpos (image ->£p, fcimage- >fpos ); 

for (i-0/i<YDIM/i++) { 

fgetpos ( image ->fp, &tmp fpos ) ; 

f read (image- >buf, sizeof "{unsigned char) , XDIMR, image- >fp) ; 
pi mage = image- >buf; 
for (j«0; j<3CDIM; j++) { 

dtemp = ( (double) *pimage - dlowest) *scale; 

ifidtemp < 0.0) * (pimage++) = 0; 

else if (dtemp > 2 55 . 0 ) * (pimage++) « 255; 

else *(piraage++) e {unsigned char) dtemp; 

f setpos (image- >fp, &tmp__fpos ); 
f write (image- >buf , sizeof (unsigned 
char) ,XDIMR, image* >fp) ; 

} } 



f setpos (image- >fp, &iraage->f pos } ,- 
f setpos (cortex. fp, &cortex.fpos ); 
for <i-0;i<YDIM/i++) { 

f read (image- >bu£, sizeof (unsigned char ), XDIMR, image- >fp) ; 
^ f write (image- >buf, sizeof (unsigned char) f XDIMR, cortex. fp) ; 

return ( 1 ) ; 

} 



int c lea r__int_image (struct int_buf *wordbuf f er) { 
long i , j ; 
int *pword ; 
f pos_t tmp_fpos ; 

f setpos (wordbuf fer- >fp, fcwordbuf f er- >fpos ); 
for (i=0;±<YDlM;i++) { 

fgetpos (wordbuf fer- >fp, &tmp_fpos ); 

pword ■= wordbuf f er - >buf ; 

f read (wordbuf fer- >bu£, sizeof (int) , XDIMR, wordbuf fer->fp) ; 
for (j-O; j<XDIM,- j++) * (pword++) « 0; 
f setpos (wordbuf fer- >fp, &tmp fpos ); 

fwrite (wordbuf fer->buf, siieo?( int) ,XDIMR, wordbuf fer->fp) ; 
return (1) ; 

} 

double find_mean_int (struct int_buf * wordbuf fer) { 
long i , j ; 
int *pword; 
double meanoO.O; 

f setpos (wordbuf fer- >fp, &wordbuf f er- >f pos ); 
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f or (i-0 ; i«eYDIM; i++> { 

pword - wordbuf f er->buf 

f read(wordbuff er->buf , eizeof (int) , XDIMR, wordbuf fer- >fp> ; 
for < j«0; j<XDIM; mean (double) * (pword++) j 

mean /= ( (double) XDIM * (double ) YD IM) ; 

return (mean) ; 

} 

int add_uchar_to_int (stmct uchar^buf * image, struct int_buf *word) { 
unsigned char *p image ; ~* 
int *pword; 
long i , j ; 
fpos — t trap_fpos; 

fsetpos (image- >fp, &image->fpos ); 
f setpos (word->xp f &word->fpos ); 
for (i=0;i<TDIM;i-*-+) { 

pword - word- >buf ; 

f getpog (word->fp, &tmp__fpos ); 

fread (word- >buf ,sizeof( int) , XDIMR, word- >fp) ; 

p image » image ->buf ; 

fread { image- >buf r e i zeof (unsigned char) , XDIMR, image- >fp) ; 
f or ( j-0/ j<XDIM/ j++) * (pword++) +- (int) * (pimage++) ; 
fsetpos (word->fp, &trnp_fpos j; 

f write (word- >buf , sizeof ( int) , XDIMR, word- >fp) ; 
return (1) ; 

} 



int add_char__to_uchar_crea t ing_uchar ( s true t char_buf *cimage, 
struct uchar_buf * image, 
struct uchar_buf *out_itnage) { 
unsigned char *p image, *pout_image ; 
char *pcimage; 
int temp; 
long i , j ; 

fsetpos (image- >fp, &image->£pos ) 

f setpos (out_image- >fp, iout_image- >f pos ); 

f setpos (cimage->fp, tcimage- >f pos ); 

f or (i-0;i<YDIM;i++) { 

pc image — cimage->buf / 

fread (cimage- >buf , si zeof (char) , XDIMR, c image ->fp) ; 
p image » image- >buf ; 

fread (image- >buf, sizeof (unsigned char) , XDIMR, image- >f p) ; 
pout_image e= out_ image- >buf ; 
f or (j=0; j<XDIM; { 

temp = (int) * (pimagef+) + (int) * (pcimage+f ) ; 

if (temp<0) temp = 0; 

else if (temp > 255) temp » 255; 

* (pout — image++) « (unsigned char) temp; 

f write (out_image->buf, sizeof (unsigned 
char) , XDIMR, out_image- >f p) ,- 

return (1) ; 

} 



int copy_int_to_int (struct int_buf +word2 , struct int_buf *word) { 
long i ; 

f setpos (word2- >fp, fcword2->fpos ); 
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f set pas (word- >fp, &word->fpoe ); 
for (i-0;i<YDIM ;£+-*■) { 

f read (word- >buf , si2eof ( int J , XDIMR, word- >f p) ; 

f write (word- >buf , sizeof (int) , XDIMR, word2 - >fp) ; 

return (l) ; 

} 



void get_snow_images (void) { 

unsigned char *psnow, *pterap; 
int number_snow__inputs ; 
int temp, *pward, *pword2, bit ,- 
long i, j ; 
double rms , dtemp ; 

live_video ( ) ; /* device specific */ 

printf ( n \n\nPlease point camera at a medium lit blank wall. n ) ; 
printf ( "\nDefocus the lens a . bit as well ") ; 

printf ( n \nlf possible, place the camera into its highest gain, 
and " ) ; 

printf ("\nput the gamma to 1.0.") ; 

printf (" Ensure that the video is not saturated ") ; 
printf ( "\nPresa any key when ready. . . 11 ) ; 

while ( ikbhitO ); 

printf ("\nNow finding difference frame rms value... ") ; 

/* subtract one image from another, find the rms difference */ 

live_video ( ) ; 

wait_vertical_blan)cs (2} 

grab_f rame ( tterap_image) ; 

live_video ( ) / 

wai t_vertical blanks <2) ; 

grab_f rame (&dTstributed_image) ; /+ use first image as buffer */ 
rms - 0.0; 

f setpos (temp^image . f p , tterap_image . fpos ) ; 

fsetpos (distributed image . fp . ^distributed image. fpos ); 

for <i = 0;i<YT>IM;i + + )7' 

ptemp = t emp_i mage .buf ; 

f read ( temp_image . buf , sizeof (unsigned 
char) , XDIMR, temp^image . fp) ; 

psnow m d±stributed_image . buf ; 

fread (distributed_image .buf , sizeof (unsigned 
char) , XDIMR, distribut ed__i mage . fp) ; 

for ( j^O; j<XDIM; j++) { 

temp = (int) *{psnow++) - (int) * (ptemp++> ; 
dtemp = (double) temp ; 
dtemp *= dtemp ; 
reus dtemp; 



rms /- ( (double) XDIM + (double) YDIM ) ; 
rms - eqrt ( rms) ; 

printf {"\n\nAn rms frame difference noise value of %lf was 
found. " , rms) ; 

printf ( ,r \nWe want at least %lf for good measure " , RMS VAJL) ; 
/* we want rms to be at least RMS_VJLL DN, so . . . *7 
if (rms > RMS_VAL ) numbe r_s no w_ inputs = i; 
else { 

dtemp = RMS__VAii / rms; 

dtemp *= dtemp; 
^ number _s no w__inpute = l + (int) dtemp; 

printf ( M \n%d images will achieve this noise 
level M , numbe r_s no w__xnputs ) ; 
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/* now create each snowy image */ 

printf ( "\nStarting to create snow pictures... \n") ; 
f setpos (snow_imagec . fp, fcsnow_images . fpos ); /* set on first 
image*/ 

for(bit = O; bit < BITS; bit++) { 

clear_int_image (ttemp_wordbuf f er) ; 
for (i=0;i<number_snow_inputs ; i+ + ) { 

live_ video () ; 

wait_vertical_blanfcs (2) ; 

graJb_ frame (&temp_image) ; 

add_uchar_to_int (&temp_image, ttemp_wordbuf f er) ; 

clear_int_image ( & temp_worc3buf f er2 ) ; 
for (i»0; i<number_snow_inpuCS ; i+*+) { 

live_video() ; 

wait_vertical_blanJcs (2) ; 

grab_f rame (&temp_ image) ; 
j add_uchar_to_int < &temp_image , &temp_wordbuf fer2 ) ; 

/* now load enow_images [bit] with the difference frame 
biased by 

123 in an unsigned char form just to keep things clean */ 
/* display it on cortex also */ 

f setpos ( temp_wordbuf f er2 . fp, &temp_wordbuf f er2 . fpos ) ; 
f setpos (terap_wordbuf f er . f p, &temp__vordbuf f er . fpos ) ; 
f setpos ( temp_ image . fp, ttemp image . fpos ) ; 
f or (i«0; i<YDiM;i++) { ^ 

pword =» temp_wordbuf f er.buf ; 

f read ( temp_wordbuf f er . buf , a izeof ( int ) , XDIMR, temp_wordbuf 



fer.fp) ; 
f f er2 . fp) 



pword2 = temp_wordbuf f er2 .buf 

f read ( terap_wordbuf f er2 . buf , s izeof (int) , XDIMR, temp_wordbu 



psnow = snow_images . buf ; 
ptemp m temp_iraage .buf ; 
for ( j-0/ j<XDIM; j+f) { 
. , ^ Mpsnow++) = * (ptetnp-n-) » (unsigned char) 

(*(pword++) - * (pword2++) ■+ 128); 

f write (snow_images.buf , sizeof (unsigned 
char) , XDIMR, s no w_ images . f p) ; 

f write (terap_image . buf , sizeof (unsigned 
char) , XDIMR, temp_ image . f p) ; 

f reeze_f rame ( ) ; 

display_uchar < &temp_ image , 0) ; /*x signifies to stretch the 
contrast*/ 

printf ( M \rDone snowy %d n ,bit); 
^ wait_vertical_blanks (30) ; 

return; 

} 



void loop_visual (void) { 

unsigned char *psnow; 

char *pcomp; 

long i , j , count «■ 0 ,- 

int ok«0, temp, bit, add__it; 
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double scale = 1.0/ RMS_VAL / 
doubl e dtemp , tmps cal e ; 
f pos_t trap_f pos ; 

/* initial rms of each, snowy image should be around 5 to 10 DN. 
let's assume it is 5, and assume further that our acceptable 
noise level of 

the full snowy composite should be approximately 1 DN, thus we 
need to 

scale theru down by approximately S*BITS as a first guess, then 
do the 

-visual loop to zoom in on final acceptable value * / 

printf ("\n\n Now calculating initial guess at amplitude... 

while ( ! ok ) { 

/* calculate anow^compoaite */ 

/* clear composite */ 

cl ear_char^image (&snow_compo3ite) ,* 

f setpos (snow^images . fp f &snow images . fpos } ,- /* set on 
first image*/ 

for (bit=0;bit<BITS;bit++) { 

j = 128 :>> bit; 

if ( our_code & j)add_it«l; 

else add_it=0,- 

f setpos (snow^composite . f p r tsnow composite . f pos ) ; 
for <i~0;i<YX>IM;i++) { " 

psnow «a snow_images . buf ; 

f read ( snow_images . buf , sizeof (unsigned 
char) , XDIMR , snow_ images . fp) ; 

fgetpos (snow_composite. fp, &tmp_fpos ); 

f read ( snow_composi te . buf , sizeof (char) , XDIMR, snow_coraposite . fp) ; 
pcomp - snow_composite.buf ; 
for (j-0 , j <XDIM; j ++) { 

dtemp - < (double) * (psnow++) -128.0) * scale; 
if (dtemp<0 , 0) { 

temp - - (int) f abs ( -dtemp + 0.5),- 

else { 

temp = (int) f abs ( dtemp +0.5); 

if(temp 127) { 
temp = 127; 

else if (temp < -12B) { 

temp = -128; 

if (add_it) { 

*(pcomp++) +» (char) temp; 

else { 

* (pcomp++) -= (char) temp; 

> } 

f setpos (snow_composite . fp, &tmp_fpos ); 
f write (snow_composite. buf , sizeof (char) , XDIMR, snow_composite . fp) ; 
^ printf ("\rDone snowy Vd tf ,bit); 

/* add snow composite to test image to form dist image */ 
add_char_to_uchar_creating_uchar ( 

&snow_composite , 

&test_image , 
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£diefcrife\ifc&d_imago) ; 

/* display botH and cue for putting scale down, up or ok */ 
i= count o 0; 

printf ( n \n Depress any key to toggle, enter to move on . . . \n 

M ) ; 

print f C"\r Distributed Image. . . ") ; 

display_uchar ( &distributed_iraage, 0) ; 
while ( getcht) 1= ' \r' ){ 
if C (couat++) % 2) { 

prititf ("\r Distributed Image... " ) ; 

display__uchar < &dis tributed — image , 0 ) ; 

else { 

printf ("\r Original Image... n ) ; 

display__uchar (&test_image , 0 ) ; 

printf { a \nScale - %lf », scale); 

printf {"\nEnter new scale, or >leC for ok... •» ) ; 
scanf < • Vlf * , & tmpscale) ; 
if (tmpscale > le6)ok«*l; 
else scale = tmpscale; 

/* distributed image now is ok; calculate actual snow_images 
used and 

store in those arrays; */ 

f aetpoa {anow_irnages .f p, tsnow_iraagea . fpos } ; /* set on first 
image*/ 

printf ( M \nNow storing snow images as used... \n" ) ; 
for <bit«0;blt<BITS;bit++) { 
for (i-0;i<YDIM;i++) { 

psnow = saow_images.buf ; 
fgetpos (snow_imagea . f p, &tmp_fpos ); 
fread (snow_image3 .buf , si2eof (unsigned 
char ) , XD I MR , snow image s . f p ) ? 

for (J-O; j <XDIM; j ++) { 

dtemp *» ( (double) *psnow -128.0) * scale; 
if (dtemp<0 . 0 ) { 

temp » - (int) f abs < -dtemp +0.5); 

else { 

temp = (int) f abs ( dtemp +0.5] ; 

*(psnow++) a (unsigned char) (temp + 128); 

f setpos (snow_images . fp, &tmp_fpos ); 
f write (snow_images.buf , sizeof (unsigned 
char) ; XDIMR, snow_images . fp) ; 

printf ( * \rDone snowy %d " ,bit); 

return ; 

> 



double find__gr ad (struct int_buf * image, int load_buf fer2) { 
int buf 1 [XDIMR) ,buf2 [XDIMR] ,buf3 I XDIMR) ; 
int +pbuf 1 , +pbuf 2 , *pbuf 3 , *p2 ; 
double total =0 . 0 , dtemp ; 
long i, j ; 
f pos_t tmp_pos ; 
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fgetpos ( image- >fp, 
fgetpos (image- >fp, 



&image - >f pos 
Sktntp_ pos ) ; 



f setpos (temp_wordbuf fer2 . fp, & t emp__wo r dbu f f er 2 . fpos > ; 

for (i»l;i< (YDXM-l) ;i++) { 

f setpos (image- >fp, stmpjjos ); 

f read (buf 1, sizeof (int) ,XDIMR, image- >fp) ? 

fgetpos (image- >fp, &tmp_ pos ); 

f read (buf 2 , sizeof (int) ,XDIMR, image- >fp) ; 

f read (buf 3 , sizeof (int) , XDIMR, image- >fp) ; 

pbuf l*=buf X; 

pbuf2=»buf 2; 

pbuf3=buf 3; 

p2 « temp_wordbuf f er2 .buf ; 

if (load_buf f er2) { 

for(j«l; j<(XDIM-l) ; { 

dtemp - (double) * (pbuf 1++) / 



dtemp 
dtemp 
dtemp 
dtemp - 
dtemp 
dtemp •+- 
dtemp + 
dtemp + 
*p2 o 
if( 

} 

else 



(double) * (pbuf 1++) 
(double) * (pbuf 1 — ) 
(double) * (pbuf 2++) 
(8.0 * (double) * (pbuf 2++) ) 
(double) * (pbuf 2--) 
(double) * (pbuf 3++) 
(double) * (pbuf 3++) 
(double) * (pbuf 3--) 
(int) dtemp; 
p2 > GRAD_THRESHOLD ) { 
(p2++) QRADJTHRESHOLD ; 



} 



if ( *p2 < -<3RAD_THRESH01*D ) { 
Mp2++) += GRAD_TKRESHOLiD ; 

else { 

*(p2++) - 0; 

} 



f write (temp^wordbuf fer2 . buf , sizeof (int) , XDIMR, temp__wordbuf f er2 . f p) 
else { 

f read (temp_wordbuf f er2 . buf, sizeof (int) , xdimr, temp_wordbuf f er2 . fp) ; 
for (j=l; j< (XDIM-l) ; j++) { 

dtemp » (double) * (pbuf 1++) ; 
dtemp (double) * (pbuf 1+-0 
dtemp += (double) * (pbuf 1--) 
dtemp += (double) * (pbuf 2+-*-) 
dtemp -a (8.0 * (double) * (pbuf2++) ) 
dtemp -t-= (double) * (pbuf 2--) 
dtemp -♦-= (double) * (pbuf 3++) 
dtemp (double) * (pbuf 3++) 

dtemp += (double) * (pbuf 3--) 

dtemp -« (double) * (p2++) ; 

dtemp *«= dtemp; 
total += dtemp ; 



return (total) ; 
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void search 1 (struct uchar_buf +suspect) { 
unsigned char *psuspact , *psnow; 
int bit, *pword,temp; 
long i , j ; 

double add^raetric , subtract_metric ; 
fpos_t tmp_fpos; 

/* this algorithm is conceptually the simplest. The idea is to 

step 

through each bit at a time and merely see if adding or 
subtracting the 

individual snowy picture minimizes some 'contrast' metric. 
This should be the most crude and inefficient, no where to go 

but 

better */ 

f setpos (snow_images . fp, &snow_images . fpos ) ; 
temp=256 / 

clear_int_image ( &temp_wordbuf f er) / 
add_uchar to^int (suspect, &ternp_wordbuf f er ) ; 

f ind_grad (&temp_wordbuf f er , 1) ; /* 1 means load temp wordbuf f er2 

*/ 

for (bit»0;bit<BITS ;bifc++) { 
/* add first */ 

f getpos (snow_iraages . f p, &trap_fpos ) ; 
f setpos (suspect- >fp, &suspect- >f pos ); 

f setpos ( t emp_wor dbuf f er . fp , £temp_w©rdbuf f er . fpos ); 
for <i«=0;i<YDIM;i++) { 

pword » terap_wordbuf f er . buf ; 
psuspect = suspect- >buf ; 
psnow = snow_images .buf ; 
f read (suspect- >buf , sizeof (unsigned 
char) , XDIMR , suspect ->fp) ; 

f read (snow_images .buf, sizeof (unsigned 
char) , XDIMR, snow images . fp) ; 

for (3»0,- j<XDIM; j++) { 

* (pword++) « (int) * (psuspect+f) + (int) * (psnow+ + ) -128 ; 

f write (temp_wordbuf f er .buf , sizeof (int) , XDIMR , temp wordbu 
ffer.fp); ~~ 

add_roetric - f ind_grad ( &temp_wordbuf f er , 0) ; 

/* then subtract */ 

f setpos ( s no w_ images . fp, ttmp_£pos ); 
f setpos (suspect- >fp, &suspect->£pos ); 

f setpos (ternp_wordbuf fer. fp, &temp_wordbuf f er . fpos ); 
for (i=0;i<YDIM;i++) { 

pword = temp_wordbuf f er .buf ; 

psuspect » suspect- >buf ; 

psnow = snow_images . buf ; 

f read (suspect- >buf , sizeof (unsigned 
char) , XDIMR. suspect ->fp) ; 

f read (snow_imagea . buf , sizeof (unsigned 
char) , XDIMR, s no w_i mages . f p) ; 

for (j = 0;jOCDIM; { 

^ * (pword+t) « (int) * (psuspect++) - (int) * £psnow++) -+12 8 ; 
f write (temp_wordbuf fer .buf , sizeof (int) , XDIMR , temp_wordbu 

ffer. fp) ; 

subtract_metric = f ind_grad ( &temp_wordbuf fer , 0) ; 

printf ("\nbit place %d: add-%le , 
sub»*fle" ,bit , add_metric, subtract_metric) ; 
temp/ =2 ; 

if (add_metric < subtract_metric) { 
printf (" bit value * O" ); 
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else { 

printf (" bit value « 1") / 
found__code +■ temp; 

} } 

printf ("\n\nYour magic number was %d" , found code) ; 
^ return ; ~ 



void search_2 (unsigned char * suspect) { 
if (suspect) ; 
return ; 



void loop__ simulation (void) { 

unsigned char *ptemp, *pdist 
int *pword, int mean , ok-O , temp ; 
long i, j ; 

double mean, scale ; 

/* grab a noisy image into one of the temp buffers */ 
printf ( " \ngrabbing noisy frame. . . \n") ; 
clear_int_image ( &terap_vordbuf £ er) ; 
for (i=0;icNOM_NOISY/i++) { 

live_video() ; 

wait_vertical_blanks <2) ; 

grab_ frame ( &temp_image) ; 

add uchar_to__int <& temp image, & temp wordbuf f er) ; 
j - (Tong ) NUM NOISY ; " ~ 
printf (»\r%ld of %ld ",i+l,j) ; 



/* find mean value of temp__wordbuf f er */ 
mean - f ind_mean_int (ttemp_wordhuf f er) ; 
in t__mean = ( int ) mean ; 

/* now we will add scaled version of this 'corruption' to our 
distributed 
image */ 
scale c= 1.0; 
while ( iok ) { 

/* add noise to dist image storing in temp_ image */ 

f setpos (distributed_image . fp, &distr ibuted_image . f pos ) ; 

f setpos ( temp_wordbuf f er . fp, &t empower dbuf f er . fpos ) ; 

f setpos (temp_image . fp, &temp image. fpos ); 

for (i=0;i<YDIM;i+^) { ~ 

pdist = distributed_image .buf ; 
pword = temp_wordbuf f er.buf ; 
ptemp = temp_iraage.bu£j 

f read (dia tr ibuted_image , buf , eizeof (unsigned 
char) , XDIMR, dis tribut ed_image . fp) ; 

fread(temp wordbuf f er . buf , s i zeo f ( int ), XDIMR , temp wordbuf 
fer.fp); ~ 

for ( j=0; j<X0IM; { 

temp m (int) *(pdist+-f) + * (pword++) - int_mean; 

if (temp<:0) temp - 0; 

else if (temp > 25 5) temp » 255; 

* (ptemp ++) = (unsigned char) temp/ 



- 71 - 



2000 09 11 13:32 



9-509795 



fwrite(temp image .buf, sizeof (unsigned 
char) , XDXMR , temp_image . rp) ; 

/* display the dist image and the corrupted image */ 
display_uchar (&temp_image, 0> ; 

/* apply new 'corrupted' image to search algorithm l for id 

value */ 

search_i (&temp_image) ; 

/* apply new 'corrupted' 'image to search algorithm 2 for id 

value */ 

/* 

search_2 (temp_image) j 
*/ 

/* prompt for upping noise content or ok */ 

OX =» 1; 

} 

return; 

} 



int initial ize_everything (void) { 
long i,j; 

unsigned char *pucbuf ; 
char *pcbuf ; 
int *pibuf; 

/* initialize cortex */ 

s t rcpy ( cortex . filename , " f : image " ) ; 

if { (cortex. fp»f open (cortex, filename, M rb") ) = =NUIjL) 
system("v f g n ); 

else f close (cortex, fp) ; 

if ( (_inp(PORT_BASE) » OxFP) ){ 

print f ("oops ") ; 

exit [0) ; 



/* open cortex for read and write */ 

if ( (cortex. fp=f open (cortex, filename, "rb+ n ) ) s=NULLJ { 
pcintf( M No good on open file joe " ) ; 
exit 10) ; 

fgetpos (cortex. fp # ficortex.fpos >; 

/* test_image; original image */ 

etrcpy (test_image . filename , "e : tst_img" ) ; 

if ( ( test_image. f p=f open (test_iraage. filename, °wb n ) ) ==NULL) { 
printf ( M No good on open file joe w ); 
exit (0) ; 

} 

pucbuf = test_image.buf; 

for (i ss 0;i<XDXM;i+-»-) * (pucbuf++) =0 ; 

f or (i=0;i<YDIM;i+-t-) { 

f write (test_image .buf , sizeof (unsigned 
char) , XDIMR, test_image. fp) ; 

f close (test_image . f p) ; 

if ( ( tea t_itnage . f p=f open ( tes t_iraage . filename f "rb+") ) ==NTXLL) { 
printf ( n No good on open file joe ••); 
exit (0) ; 
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fgetpos (test_image .fp, &test_image.fpos ); 

/* snow_composite; ultimate image added to original image */ 
etrcpy ( Gnow_composite . filename , " e : snw_cmp " ) ; 

if ( (snow_coraposite . fp=f open (snow_compo3ite . filename, "wb") ) ~=NULL) ( 
printf (" No good on open, file ioe " ) 
exit(O); 

pcbuf « snow_composite . buf ; 

for (i=0;i<XDlM;i++) * (pcbuf ++) =0; 

for(lsO;l<YX>XM;i++) { 

fwrite (snow_composite.buf , sizeof (char) , XDIMR , enow_composite . f p) ; 

f close (snov_coraposite . fp) ; 

if ( (snow_composite. f p=f open (snow_ compos ice. filename, "rb+") ) ^NULL) I 
princff" No good on open file joe ") ; 
exit(0); 

fgetpos (snow_compoaite . f p, & s no w__ composite . fpos ) ; 

/* distributed__image; test_img plus snow composite */ 
strcpy(distributed_ image, filename, "e : dst_img") ; 

if ( <distributed__image . fp=f open <distributed__image . filename, "wb" ) ) ==NU 

LiIj) { 

printf ( " No good on open file ioe 11 ) - 
^ exlt(O); 

pucbuf m distributed__image_buf ; 
for (i=0 ; i<XDIM; i++) * (pucbuf++) =0 ; 
for<i=0;icYDIM;i++) { 

fwrite (distributed_image . buf , sizeof (unsigned 
char) , XDIMR , distributed_image . f p) ; 

f close (distributed_image . fp) ; 

if ( (distributed_image . f p-f open (dis tributed^itnage . filename , "rb+") ) «n 

UXjLiJ { 

printf ( " No good on open file ioe M ) ; 
exit{0); 

fgetpos (dis tributed_image . f p, &distributed_image . fpos ) ; 

/* temp_image; buffer if needed */ 

st rcpy (temp_image , filename ,» e : temp_img ") ; 

if ( Ctemp_image . f p=f open <temp_xmage. filename, "wb") )==NDIiL) { 

printf (" No good on open file joe ■); 
^ exit ( 0 ) ; 

pucbuf = temp_ image . buf ; 

for ( i*= 0 ; i<XDIM ; i++ ) * (pucbuf ++) =0 ; 

for (i=0;i<YDIM;i++) { 

fwrite ( temp_image .buf, sizeof (unsigned 
char) , XDIMR # temp_iraage . f p) 

f close (temp^iraage. fp) ; 

if ( (temp_image . fp-f open (temp_image . filename, "rb+") ) «=NULL) { 
printf ( " No good on open file ioe ") ; 
exit (0) ; 

} 

fgetpos ( temp_image . f p , &temp_ image . fpos ) ; 
/* temp_wordbuf fer; 16 bit image buffer for averaging */ 
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Btrepy <temp_wordbuf f er. filename, "e: temp_wrd r ) ; 

if ( ( temp_vordbuffer . f p»f open (temp_wordbuf f ear. filename, " wh " ) > ==NULU) { 
printf ( " No good on open file joe " ) ; 
exit (0) ; 

pibuf = temp_wordbuffer . truf ; 
for <±«0;i<XDIM;i++) * (pibuf ++) »0; 
for <i=0;i<YDIM;i++) { 

fwrite (temp_wordbuf f er.buf , sizeof (int) , XDIMR, terap_wo:rdbuf f er . f p) ; 

f close (temp^wordbuf fer . fp) ; 

if ( ( temp wordbuf f er . fp«f open (temp wordbuf f er . filename , " rb+ ■ ) ) «— NUIjL) 
{ ~~ 

printf (" No good on open file joe 11 ) ; 

exit <0) ; 

fgetpos (temp_ wordbuf f er . fp, ttemp_wordbuf fer . f pas ); 

/* temp__wordbu£ f er2 ; /* 16 bit image buffer for averaging */ 
strcpy ( temp_wordbuf f er2 . filename , M e : tmp_wrd2 " ) ; 

if < (temp_wordbuf fer2 . fp=f open (temp wordbuf fer 2 . filename , "wb" ) ) —NULL. 
){ 

printf ( " No good on open file joe " ) ; 
exit (0) ; 

pibuf * temp_wordbuf f er2 . truf ; 
for (i=0;i<XDIM;i++) * (pibuf ++) =0; 
for (i=0;i<YDIM;i++) { 

fwrite ( t e mp_wordbuf f e r 2 .buf , sizeof (int) , XDIMR, temp_wordbuf f er2 . f p) 

f close (temp_wordbuf f er2 . fp) ; 

if ( (temp wordbuf fer2 . fp- fop en ( temp wordbuf £er2 - filename, "rb+") ) »«NUL 
D { 

printf!" No good on open file joe ") ; 
exit (0) ; 

} 

f getpos (temp_wordbuf £ er2 . fp , &temp_wordbuf f er2 ■ fpos ) ; 

/* snow^images ; BITS number of constituent snowy pictures +/ 
strcpy (gnow^imagea . filename, " snw_imgs n ) ; 

if ( (snow_ images . f p=f open (snow__ images . filename , "wb") ) x==NUX»L) { 
printf (" No good on open file joe 
exit (0) ; 

} 

pucbuf « enow_images . buf ; 

f or (i=0;i<XDIM;i++) * (pucbuf ++) »0; 

f or <j = 0; j<BITS; j++) { 

f or (i=0;i<YDIM,-i + + > { 

fwrite (snow_imagea .buf , sizeof (unsigned 
char) , XDIMR, snow_images . f p) ; 

f close (gnow_i mages . fp) ; 

if ( ( snow_images . f p~f open (snow_images . filename, " rb+ u ) ) «==NULi!L) { 
printf { D No good on open file joe n ) ; 
exit <0) ; 

f getpos (snow_images , fp t tsnow_images . fpos ); 
return (1) ; 
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int close_everything (void) ( 

fcloee (teet^image.fp) / 

f close (snow^coraposite . f p) ; 

f close (distributed_image . f p) ; 

f close ( t emp^image - fp) ; 

f close (ternp_wordbu£f er. fp) ; 

f close (cemp_wordbuf f er2 . fp) ; 

f close (snow_images . fp) ; 

return (1) ; 

} 



main () { 

int i, j; 



V 



printf ( "\nlnitializing. . .\n\n") ; 

initialize_everything 0 ; /* device specific and global malices 
live video { ) ; 



/* prompt for which of the three demos to perform * / 
while ( demo <. 1 j | demo > NUM_DEMOS } { 

print f ( "Which demo do you want to run?\n\n") ; 
print f ("1: Digital Imagery and Very High End Photography 
Simulation\n") ; 

printf ("2: Pre-exposed Print Paper and other Dupping\n") ; 
printf ("3: Pre-exposed Original Film (i.e. In- Camera) \n") ; 
printf ( "\nEnter number and return: w ) ; 
scanf ( " td" , &demo) ; 

if (demo c 1 | | demo > NUM_DEMOS3 { 
printf ( M \n eh eh •■) ; 

} } 

/* acquire test image */ 

printf ("\nPress any key after your test scene is ready... ") ; 
getch ( ) ; 

grab__f rame {&tast_image) ; /*grab_f rame takes care of device 
specific stuff*/ 

/* prompt for id number, 0 through 255 */ 
printf ( "\nEnter any number between 0 and 255. \n n ); 
printf ("This will be the unique magic cade placed into the 
image : " ) ; 

scanf ( w *d n , &our_code) ; 

while (our_code<l | f our_code>25S) { 

printf ( * Between 0 and 255 please ") ; 

scanf ( Hd" , Sour code); 

} 

/* feed back the binary code which will be embedded in the image 

printf ("\nThe binary sequence "); 
for (i=0;i<BITS;i++) { 
j o 128 >> i; 

if( our^code & j ) printf ( n 1 ") ; 
else printf ("0" ) ; 

printf ( n (%d) will ba embedded on the image\n" , our_code) ; 

/* now generate the individual snow images */ 
get_snow_images () ; 
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loop_visual {) ; /* this gives visual feedback on 'tolerable' 
noise level */ 

print f ( "\nWe' re now to the simulated suspect... \n n ) / 
loop_simulation < ) ; 

close_everything ( ) ; 
return (0) ; 
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